A reinterpretation of the metamaterial perfect absorber 
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We analytically prove that the appearance of two almost, but not exactly out-of-phase currents inside a 
metamaterial is necessary for a perfect absorber. We further show that although evanescent waves do not 
contribute to far fields scattered by the perfect absorber, they however must consume the electromagnetic 
energy significantly. © 2012 Optical Society of America 
OCIS codes: 160.3918, 250.5403 



Iii 2008, Landy et al. experimentally demonstrated 
that a metamaterial consisting of three thin metal- 
dielectric-metal layers can almost entirely absorb an elec- 
tromagnetic (EM) wave at particular frequencies, which 
they called the metamaterial perfect absorber (MPA) 
[1]. The same group [1-4] also addressed the underly- 
ing mechanisms of the MPA by using the popular ef- 
fective bulk permittivities and permeabilities employed 
in the metamaterial community [5]. It was found that 
the retrieved effective permeability can be described by 
a Lorentzian function around the designed frequency, in- 
dicating the appearance of a magnetic resonance [4] . This 
observation was further supported by full-wave simula- 
tions which showed the currents in the two metallic layers 
form a circulating current loop [3] . 

Recently, Chen et al. proposed a different interpreta- 
tion [6,7]. They noticed that there are two interfaces in- 
side a three-layer metamaterial, and EM waves are scat- 
tered and reflected by these discontinuities. It is therefore 
natural to treat these layers separately and further take 
into account multiple reflections inside the metamate- 
rial [8,9] . By treating each metallic layer as a metasurface 
with effective surface electric and magnetic susceptibili- 
ties [10], their model can provide almost identical spectra 
as that obtained by full-wave simulations. Furthermore, 
the currents in the two metallic layers were observed to 
have almost equal amplitude and to be nearly out-of- 
phase [7]. 

In this paper, we reinvestigate the underlying physics 
of MPAs by studying their polarization currents. We find 
that MPAs possess not only a magnetic dipole but also 
an electric dipole and quadrupole as well. These mul- 
tipoles are destructive along the reflected direction and 
constructive along the transmitted direction but out-of- 
phase with the external incident wave. We further show 
that non-zero-order Bragg waves are crucial for the EM 
energy consumption. 

Let us consider a metamaterial membrane, surrounded 
by vacuum, which is periodic in the xy plane. For sim- 
plicity, we assume its metaatoms are arranged in a rect- 
angle lattice with primitive lattice vectors d x e x and 
d y e y . Here d x and d y are the corresponding lattice con- 
stants. Wc further assume that the external illumina- 



tion is a plane wave propagating along the z direction 
with a wave vector k z . Under this external excitation po- 
larization currents, J(r, oj) = — iwe [e(r,w) — l]E(r,w), 
will appear inside the metamaterial, where e(r,w) rep- 
resents the permittivity of the metamaterial. Using the 
free space Green's function, we can express the scattered 
field in terms of the current 

E(r>) = ]T ~T^° e ' k "'" r / drJ„ m , ± (r, W )e- lk -- r , 

where the integration is performed over one unit cell, 
Zq = sj Ho/ e o is the free-space intrinsic impedance, A is 
the excitation wavelength, and fco = 2tt/X is the free- 
space wave number. The wave vector k mn is defined as 

kmn = k|| gmn i ^rnn^z with K<rnn = \J |^mn,|| 

and g m „ = {2nm/d x )e x + (2Trn/d y )e y being the re- 
ciprocal wave vectors, and k|i being the projection of 
the incident wave vector onto the xy plane. Here the 
positive sign corresponds to forward scattering (prop- 
agating along the positive z direction), and the nega- 
tive sign corresponds to backward scattering (propagat- 
ing along the negative z direction). Moreover, J m „,i = 
J — k mn (k mn -J)/fcQ. Consequently all orders of scattered 
waves are transverse. Under certain circumstances, only 
the zero order waves can survive in the far zone. The 
forward wave is then given by 

Eu») = J^^ V ' \ rfrJ 00 , ± (r, W ) e - kI -. 

The backward wave bears a similar expression except 
that k* is replaced with k|| — k\e z . This equation sug- 
gests that the tangential component of the incident wave 
vector is conserved. 

Wc now apply these equations to the metamaterial 
perfect absorber. The most studied MPA is the one pro- 
posed in Ref. [4] (depicted schematically in Fig. 1), which 
consists of three thin layers. The first layer is a perforated 
metallic membrane which is referred to as the cross layer. 
The second layer, called the spacer layer, is filled with 
a homogeneous dielectric medium with weak absorption, 
and the last layer is generally a metallic ground plane. 
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Furthermore, the meta-atoms of the cross layer are ar- 
ranged periodically in a square lattice. The whole struc- 
ture possesses both x and y mirror symmetries. Under 



normal incidence W 



„ik z' 



e x , we know by symmetry 



that E x is an even function of x, while E y and E z are 
odd functions of x. Furthermore, the lattice constant d 
is far smaller than the incident wavelength, so that only 
the zero-order waves survive in the far zone. The forward 
and backward scattered fields are hence given by 



h/2 



g{z,uj)e- lk0Z dz 



-h/2 
h/2 



h/2 



g{z,uj)e ik ° z dz. 



Here h is the thickness of the MPA, and the integral 
function 



g(z,u) 



1 

d 2 



d/2 r-d/2 



-d/2 J -d/2 



J x (x,y,z,u)dxdy (1) 



stands for the current density at a specific z plane. Note 
that J v does not contribute because of its symmetry. As 
a direct result, the reflected wave is E b , and the trans- 
mitted wave is 

A metamaterial with unity absorption neither reflects 
nor transmits the incident EM wave, which is equivalent 
to requiring 



h/2 

g(z)e ik ° z dz = 0, 

-h/2 

h/2 r, 

g(z)e- lkoZ dz 

-h/2 



Z ' 



(2) 
(3) 



To simplify the notation, we dropped the lu dependence 
of g here and in the following. The first equation imposes 
no reflection, and the second one leads to zero transmis- 
sion. Because the thickness h is much smaller than A, 
we can employ multipolar analysis by expanding e tk ° z 
as 1 + ikoz + • ■ • . The leading order gives the electric 
dipole and the first order of k$z gives the combination 
of magnetic dipole and electric quadrupole. Supposing 
these three multipoles dominate the above two integra- 
tions, we arrive at 



,fc/2 i rh/2 

/ 9R(z)dz=—, / gi(z)dz = 0, 

J-h/2 Z J-h/2 



h/2 
-h/2 



(4) 



as well as 

,h/2 



rh/2 1 
g R {z)zdz = 0, k g I (z)zdz = —, (5) 

h/2 J-h/2 Z 



where gn and gj are the real and imaginary parts of 
g, respectively. Eqs. (4) immediately suggest that there 
exists an electric dipole moment with a value of 1/Zq, 
and Eqs. (5) imply that the MPA possesses a magnetic 
dipole as well as an electric quadrupole, and their sum 
equals i/Zo (see Eq. (9.31) of Ref. [11]). Moreover, as 
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Fig. 1. Reflection and absorption spectra of the meta- 
material perfect absorber. The inset shows the geome- 
try of the metamaterial. The metallic cross consists of 
two 0.4 x 1.7 metallic bars. The thicknesses of the cross, 
spacer and ground layer are 0.1, 0.09 and 0.2, respec- 
tively. The total thickness h equals 0.39, and the lattice 
constant d = 2. All dimensions arc in micrometers. 



suggested by Eqs. (2) and (3), these three multipoles are 
destructive along the reflected direction, which leads to 
zero reflectance. They, however, are constructive along 
the transmitted direction but out-of-phase with the ex- 
ternal incidence so that the transmittancc is also zero. 

It is important to notice that the dielectric layer gen- 
erally has a permittivity ta which is very different from 
the permittivity e m of the metallic layers. For instance, 
the ratio (e m — l)/(ed — 1) of the MPA studied below is 



about 1428e 4 ' 



/hen A = 5.93/im. Since E x is continu- 



ous crossing the spacer-ground interface, the polarization 
current J x is strongly concentrated inside the metallic 
layers. We therefore can neglect the polarization current 
within the dielectric layer. Denoting the total current of 
the cross and ground layer as G c and G g , respectively, 
Eqs. (4) suggest 



Z Q x Rc(G c + G g ) = l, Im(G c + G g 



0. 



(6) 



In addition, Eqs. (5) imply |<jj| is much bigger than \gn\ 
because k \z\ <C 1. We therefore expect that G c and G g 
are nearly purely imaginary. In other words, the currents 
in the two metallic layers are almost, but not exactly, 
out of phase. This fact, here proved analytically, was 
reported in earlier numerical simulations [1-4,12]. 

To support the statements above, let us consider a 
MPA whose geometrical and optical parameters are al- 
most identical to the ones used in Ref. [4] ; our geometri- 
cal parameters are specified in Fig. 1. The permittiv- 
ity of the dielectric is described by a Lorentz model, 
e d (uj) = too [1+lj%/(uj% - lj 2 - iwy)] with = 2.24, 
uj p = 93.77 THz, 7 = 173.73 THz and lo = 3.1 
THz [13]. The permittivity of the metal is described by 
a Drude model, e m {u>) = 1 — ui 2 m /(uj 2 + iujj m ), with 
u pm = 1.37 x 10 4 THz and j m = 40.8 THz. Using 
a finite-difference time-domain method [14], where the 
size of spatial grid cell is fixed at 5 nm, we calculate the 
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Fig. 2. (a) Amplitude and phase of the function g(z), and 
(b) the function rj{z), when A = 5.93/mi. The dielectric 
layer is highlighted. 



tribute significantly to the energy dissipation. Indeed, 
we find that the ratio of \E Z /E\ inside the spacer layer 
ranges from 0.84 to nearly 1.0. This property may have 
many important applications. For example, it may sig- 
nificantly improve optical coupling in quantum well in- 
frared photodctectors [15]. By integrating ry in each sub- 
layer, we further find that the relative contributions from 
the cross, spacer and ground layer are 28.5%, 56.7% 
and 14.8%, respectively. Therefore, the dielectric and the 
metal consume EM energy almost equally. 

To sum up, we studied the metamaterial perfect ab- 
sorber in terms of its polarization current. It is analyti- 
cally proved that there always exists one circulating cur- 
rent loop, together with an electric dipole as well as an 
electric quadrupolc. It is the interference among them 
and the external incident wave that results in unity ab- 
sorption. We further showed that longitudinal fields con- 
tribute significantly to the dissipation of the electromag- 
netic energy within the perfect absorber. 

We acknowledge support from the LANL LDRD pro- 
gram. 



linear spectra and plot the results in Fig. 1. Around a 
wavelength of 5.93 /xm, the absorption is found to be 
nearly 100%. We further calculate the current function 
g(z) at this wavelength and plot the result in Fig. 2(a). 
To check our numerical results, we computed the inte- 
grations in Eqs. (4) and (5), and verified that they are 
nearly identically satisfied. As discussed above, we find 
that the polarization current is strongly localized inside 
the two metallic layers, and a phase jump of 0.947T ap- 
pears at the spacer-ground interface. Furthermore, we 
obtain G C Z Q = (-0.15-6.77i) and G g Z = (l.l+6.78i), 
in perfect agreement with the discussions above. 

Next we consider the dissipation of the electromag- 
netic energy within the MPA. According to Poynting's 
theorem, the energy absorption can be described by J-E 
[11]. Consequently we define 

J Jlm{e)\E(x,y,z)\ 2 dxdy 
fffIm{e)\E(x,y,z)\ 2 dxdydz [ ' 

to measure the relative contribution from a specific z 
plane. Here e = e(r, uj) describes the permittivity of the 
whole MPA structure. The result is plotted in Fig. 2(b). 
It is surprising to observe that only a feeble discontinu- 
ity appears around the spacer-ground interface. Because 
of the structural periodicity, we can express the electric 
field as a superposition of Bragg waves. Under normal 
incidence, the zero-order mode, with m — and n = 0, 
represents a plane wave propagating along the z direc- 
tion. Since the longitudinal z component is absent, this 
mode must be continuous across the spacer-ground in- 
terface. We therefore expect that rj of the spacer layer 
is much smaller than that of the ground layer because 
Im(e m )/Im(ed) ~ 3000, which obviously disagrees with 
our numerical result. This contradiction suggests that 
higher-order Bragg waves or the longitudinal fields con- 
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